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Abstract. In this paper we study the large - Z behaviour of the ground state 
energy of atoms with electrons having relativistic kinetic energy sjp^c^ + m^c^ — 
mc^ . We prove that to leading order in Z the energy is the same as in the 
non-relativistic case, given by (non-relativistic) Thomas-Fermi theory. For the 
problem to make sense, we keep the product Za fixed (here a is Sommerfeld's 
fine structure constant), and smaller than, or equal to, 2/7r, which means that 
as Z tends to infinity, a tends to zero. 



1. Introduction and results 

As a model for a relativistic atom with nuclear charge Z and N electrons, we 
consider the operator 

Z ^ 1 



— . l<i<j<Af ^■'l 

Here, Xi G is the coordinate of the i'th electron, is the Laplacian with respect 
to Xi, and a is Sommerfeld's fine structure constant (the physical value of a is 
approximately 1/137.037). This is the expression one obtains using ^'p^c?- + m^c**— 
mt? for the kinetic energy of the electrons (and making the substitution p —iKV), 
measuring energies (Hrei) in units of Rydberg, and lengths (the Xi 's) in units of the 
Bohr radius. 

This model has been much studied over the past thirty years. Stability in the 
case N = 1 was proved independently by Herbst ^ and Weder The 'Stability 
of Matter' for the model was first proved by Conlon |2|, later by Fefferman and 
de la Llave 0, and also by Lieb and Yau jJHI; see the latter for an overview. A 
non-exhaustive list of other works on this model is 12UL I19L 1181 ^ . 

It is well-known that the operator Hrei is bounded from below on C§°(M.^^'^) if, 
and only if, Za < ^. Only in this case is the atom stable; and we define the 
operator H^ei as a self- adjoint, unbounded operator by Friedrichs-extending this 
semi-bounded operator. To study the energy of large atoms, one would normally 
then consider the limit as Z — > oo of the infimum of the spectrum of this operator. 
However, due to the upper bound on Z resulting from the restriction Za < — , this 
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is not possible here. To overcome this problem, we consider 

where S = Za is held fixed. This ensures that as Z — s- oo, and therefore a — > 0, 
the operator Hrei remains well-defined — as long as < <5 < f • Also, we shall keep 
A = N/Z fixed. The energy of the atom is then defined as 

En{Z,5) := MaHpiHrei) , 

where the spectrum of Hrei is calculated on Hp ~ /\^ L'^ {M.^ , C^) , the Fcrmionic 
Hilbert space, describing N Fermions, each with q possible spin states. We will take 
q = 2 from now on (but this is no restriction). We note that since (the extension of) 
Hrei is self-adjoint and bounded from below, we have the Rayleigh-Ritz principle: 
If C is a form core for the corresponding quadratic form, then 

Wee 1 11^11=1}'^ 

Our main result is the following: 

Theorem 1.1. Let S G (0, 2/7r] and X > be fixed and let Hrei and Em=\z{Zt5) 
he as above. Then 

ExziZ,5) = -CtfWZ^/^ + o{Z^/^) , Z^oo, (1.1) 

where —Ctf{X)Z'^ is the (non-relativistic) Thomas-Fermi energy of the atom. 

This shows that, to leading order, the ground-state energy of a relativistic atom is 
given by the (non-relativistic) semi-classical Thomas-Fermi energy approximation, 
as it is for the non-relativistic atom (note that the case ^ = ^ is included). (In 
the non-relativistic case this was first proved by Lieb and Simon (13j : see also Lieb 
|l()j.) This expresses the fact that for large atoms the majority of the electrons are 
non-relativistic. 

The second term in the expansion (|l.l|l will be studied in a forthcoming pa- 
per [nj. 

The proof of Theorem ll.ll will be by finding upper and lower bounds on Exz {Z,5). 
Note that the relativistic kinetic energy is always lower than the non-relativistic 
one: 

2 

\A^V+^-a"' = a"'(v/TTKF-i) <y. (1.2) 

(Note: since we will later make Taylor expansions of the square root in the rela- 
tivistic kinetic energy, we will have to insist on the non-relativistic kinetic energy 
being —A/2). This means that all upper bounds derived earlier EI| for the 
non-relativistic operator 

will also be upper bounds for Hrei] in particular, to prove Theorem 11.11 we need 
only derive a lower bound. 
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2. Organisation of the paper 

We start in SectionOlby reducing the A^-body operator Hrei to a one-particle one; 
having done that, we only need to consider wave functions given as Slater-deter- 
minants when trying to minimise the energy. To proceed, we need to localise the 
kinetic energy. To do so, we use (in Section^ an analogue of the IMS Localisation 
Formula for the Schrodinger operator, see P-27]. This formula has already been 
developed by Lieb and Yau in ^B] for both the operator V— A -I- a~'^ — and the 
hyper-relativistic kinetic energy \p\. This is essentially done by finding the integral 
kernels of these operators. For A + a^^~ a^^, this involves the modified Bessel 
function K2 , and the derivation of the formula and of needed properties of K2 are 
carried out in Appendix 1X1 The localisation error, given by a bounded operator 
L^") expressed as an integral operator involving K2, is then estimated (in Section|SJ. 
Estimating the error is rather technical (calculative) and involves localisation of the 
operator and the above mentioned properties of K2. Some of the localised terms 
are estimated with the localised energy itself (Sections and [T]). 

Coming to the localised energy, we have to estimate the kinetic energy close to 
the nucleus. Since this is the high-energy region, this is where the electrons are 
relativistic, and so this term should be of lower order, since, to leading order, there 
should be no relativistic contribution to the energy. As the relativistic kinetic energy 
is asymptotically linear in p in the high-energy region — as opposed to the classical 
one which is quadratic — the singularity in the potential causes substantially more 
trouble. This problem is solved (in Section (SJ by a clever choice of parameters in 
an estimate by Lieb and Yau in on the sum of the eigenvalues of the energy in 
a ball around the nucleus. This also determines the scale on which one can localise 
close to the nucleus. A part of two of the localised terms of the operator L*^"' is 
estimated along with this term. 

In the outer region, one uses (in Section (S)) essentially the same idea as Lieb did 
in the classical case, see ^P, to re-find the desired phase space integral, which is 
to give the semi-classical Thomas-Fermi energy. This involves introducing coherent 
states and estimating the error by doing so. The formulae for the relativistic case 
were developed in |15|. but the error obtained there is too rough for our purposes. 
We therefore develop (in Appendix^) a better estimate by a more careful analysis. 
In order to make all this work, one need the coherent state to be supported further 
out than the initial cut-off around the nucleus. To get this, an intermediary zone 
is introduced (also in Section by an additional cut-off. The energy in this shell 
is estimated (in Section [TJ by a generalised version of the Lieb-Thirring inequality, 
proved by Daubechies in Also the other part of the previously mentioned two 
terms of the localised operator L^") is estimated in this way. 

Finally we relate (in Section IS)) the energy in the outer region to the Thomas- 
Fermi energy from the classical (that is, the Schrodinger) case. In this region, the 
kinetic energy is small, and using the specific scaling property of Thomas-Fermi 
theory allows one to make the change from the relativistic energy y/^a^-^A + or^ — 
a'^ to the non-relativistic one, —A/2, getting errors of the desired order. 
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3. Reduction to a one- particle problem 
We will use the notation 



N 



H = aHrei = E { V-A, + - - ^} + (3.1) 

Recall that d = Za is fixed and that the ground state energy of Hrei is to be proven 
to be of leading order Z"^/^. Since we wish to consider a as the free parameter, the 
relevant order of all error terms will be o{a~^^^). Also, we will denote the operator 
V— A + by \/p'^ + a~^, and so T{p) = sjjp- + a"^ _ q.-i be the kinetic 
energy. 

Wc start by reducing the problem from an A'^-particlc problem to a onc-particlc 
one. This is done by using an inequality on the electron-electron interaction 
X^i<i l^i ~ ^jl "^' which will reduce this to a one-particle potential. 

Choose a spherically symmetric function g € C^(M^), non-negative, supported 
in the unit ball i?(0, 1) of M^, and such that Jg{x)'^ cPx — 1. Let (j){x) = g{x)'^ and 
let for a > (a to be chosen later), (j)a{x) = a~^<f){x/a), so that J<f)a{x)d^x = 1. 
Then for all p : ^ M we have: 

y , ' ,> y //^^fc^iM4^^rf3^d3 

l<i<j<N ^ ' ^1 l<i<j<N-'-' I 

^^1^J'' 2 JJ \x-y\ 

'til """^.f/'' A ^1/ - ^ // A 



= 1 



1 



<?io(a; - Xi) - p{x)j Y.j My - Xj) - p{y) 



2JJ \x-y\ 



d xdy 



N 



F-y| 277 \x-y\ 

In the last inequality we used that \x — y\~^ is of positive type (a positive kernel) 

since 

[[lMMMy = A.nl^d^. 
J J \x-y\ J \pr 

The constant c(^) is independent of a: 

1 ff Hx)<P{y) ,3 .3 /M^^S 

Noting that (using the spherical symmetry of 0o) 
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we get the operator inequality (see for H): 

N ^ 

-ff > 51 {yP? ^ + "P* * k^r'} (3.2) 

^JJ \x~y\ 

Having reduced the iV-body operator H to a one-body operator, we only need 
to consider Slater-determinants when trying to minimise the energy. That is, when 
considering {ijj,Hip) we need only consider those ip CzTip which are given by 

tp{xi, xn) = -^^^ det(mi(a;j)), 

where rrii e L^(R'^), i — 1, . . . ,N, are orthonormal. Note also that since C5^(]R"^) is 
a core for the operator -y/p^ ^ q,-2 _ _ (5/|a;|,(5 g [0,2/t:] (see Herbst we 
need only consider m^'s in this space. Then, as soon as /i is a one-particle operator 
acting on L^(M'^), we have that 

JV N 

i=l i=l 

Here, ( , ) and ( , ) denote inner products in L^(M'^^), respectively L^(M'^), both 
linear in the second variable, and hi is h acting on the variable Xi of tp. Also, we 
will use II • lip for the norm in LP(M'^). 



4. Localisation of the kinetic energy 

In order to treat the one-body operator in 13.2(1 and in particular the singularity 
in the Coulomb-potential — which causes considerably more trouble than in the non- 
relativistic case — we introduce, following Lieb and Yau [TB] , a partition of unity (see 
also Cycon, Froese, Kirsch and Simon [3 Definition 3.1]). For some (5 G (0, i), let 
6*1 and O2 be monotone positive smooth functions on M+, < 0^ < 1, such that 

^i(^)-\o ife>i + /? ' ^2(o-|i ife>i + /3 ' 

and such that 9i{£,f + e2{£.f = 1 for aU i e K+. Now define, with 8/9 < r < 1 and 
1/3 < t < 2/3 (these choices of parameters are governed by the later analysis), the 
following three functions, which (for a sufficiently small) is a partition of unity in 

Then, at least for a sufficiently small, we have the picture in Figure ^ 
According to Lieb and Yau | 16i Theorem 9; a^^ corresponds to m] we have for 

/ e C^(R3) that 

3 

(/, Vp^+a-2/) = X,V^+^Xjf) ~ (/, i^"V) (4.2) 
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Figure 1 . The partition of unity. 
where L^"^ is a bounded operator on _L^(M^), given by the kernel 

= ^ — ^3^2 — lJxA^)-xAy)) ■ 

Here K2 is a modified Bessel-function, defined on (0,cxd) by 

K2{t) ^7: 2;e-^*(^+^"') dx. 
2 Jo 

For completeness, we derive this in Appendix El 

Using this we find, with T{p) = + q!"^ — V{x) — ?>l\x\ and V' a Slater- 
determinant as mentioned in the previous section, that 

N 

{T(pO - V{X.) + «p * </)a * \x^\~^)i,) (4.3) 

i=l 

= ^(m„ {T(p) - t/(x) + a p * 0a * |a;r^}m,) 
1=1 

3 AT JV 

= ^^(m,,Xj{r(p) - T/(x) * |x|"^}xjm,) - ^(mj,L(")mO, 

j=i 1=1 1=1 

since X)j=i Xj(2^)^ = 1 fo^' all x G (and a sufficiently small). 

5. The localisation error 

We now estimate the error introduced by the localisation of the kinetic energy 
carried out in the last section. This error is given by a bounded operator L*^"\ 

\x,y) =\^L) \x,y) , L)'(x,y) = ^ ^-^^ ^xA^) - xM) ■ 

As noted above, this expression is derived in Appendix El We shall start by local- 
ising this operator, thereby splitting it it into twelve terms (!) which we will then 
treat individually. These terms are going to fall into groups though, and the terms 
in each of these will be estimated together by different means. Two of the terms 
will be estimated in later sections, together with the energies near the nucleus and 
in the intermediary zone, related to respectively xi and X2- 

In this section, the scale of the inner cut-off will be called that is, I — a^, 
8/9 < r < 1. Let X- be the characteristic function of the ball B(0,2l) in M.^ and 
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X+ that for the complement of this baU. Then each j ~ 1,2,3, sphts into four 
terms: 

+ X- {x)l\°'^ (x, y)x+ (y) + X- {x)lY^ [x, y)x~ {y). 
The foUowing lemma will eventually take care of six of these twelve terms: 
Lemma 5.1. Let I — a'' , 8/9 < r < 1 and assume that, with 7 = 1 — ^ > 0, 

|x| > and \y\ < hi on supp x+{x)L^j^\x,y)x-{y)- 
Then, for f e L^{R^), 

|(/,x+i5"^x-/)l <pHII/IIL 

where p{a) = o{e~'^'^" ) as a Q for all e such that < e < 7. In particular, 
p[a) = o(a") as a — > for all n G N. 

Remark 5.2. Note that the result with x and y interchanged also holds. 

Proof. By assumption we have that 

|x-y|>7|a;| on supp x+Lj"^X-- 

Since both \x\~'^ and i4r2(a~^|a;|) are decreasing in |x| (the last is obvious from the 
definition of K2), and since (xj{x) — XjivYf ^ 1, we get that pointwise, 

X^{x)Lf{x,y)x-{y) < 



-(") 
-(") 



on supp x+LjX-. Therefore 



\{f,X+L)"'x-f)\ (5.1) 



< l^J \fiy)\x-{y)d'yj (^^^ / \fix)\x+ix) 

We estimate both of these terms using the Cauchy-Schwartz inequality. For the 
first we get 

" |/(y)| x-{y)d'y < II/II2 IIX-II2 = Cp/'Wfh, (5.2) 



and for the second 



(5.3) 

Using the estimate (jA.7|) in Appendix 1X1 on K2, we get the estimate 



, K2{a-'Mx\) \\3 
X+{x) r-f^ d-'x 



< 47r 



16 7re-2" ^1^1 
]^ 2a-'^j\x\ 



(1 + {2a-^j\x\)-^ + {2a-^j\x\y^y\x\^ d\x\ 



f°° 11 
1287r2a"^7 / t'^e'^l + - + —f dt, 
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where the last equahty follows by the change of variables t — 2ja ^\x\. Dominating 
e^* in the integrand by e"*''''" and working out the resulting integral, we arrive 
at (using (|5.1|) . H5.2|) . and (|5.3|) : recall that I = a"^) 

\{Lx+Lfx-f)\ < C ||/||2a(3-5)/2^-2..'-| , 1^/2 

where 

{ . . . }'^' = {i(47rv(i-^^) + ^(47)- v(i-) 

Now, since 8/9 < r < 1, this term tends to zero as a tends to zero. Also 

^(3r-5)/2g-2^a'-i ^^(g-2.a-i) , « 0, 

for all e satisfying < e < 7. This proves the lemma. □ 

We now return to investigating the above mentioned twelve terms. Firstly, note 
that two of these terms are actually zero: 

X+{x)L'f\x,y)x+{y) = 

X-{x)L''^\x,y)x-{y) = 

as is easily seen by looking at the supports of x+i X-i Xi, and xs- Next, we note 
that the following three terms fulfill the conditions in Lemma and therefore are 
o(a"),a ^ 0, for all n G N: 

X+{x)L["\x,y)x-iy) + for > 21 and \y\ < [1 + 13)1 
X+{x)L^^\x,y)x-{y) + for \x\ > (1 - /?)a* and \y\ < 21 

X+{x)L["\x,y)x-{y) + for > (1 - and \y\ < 21 

and for \x\ e [21, (1 - ^)a*] and |y| < (1 + (3)1. 

This is due to the fact that for a small enough, a* > a'', since t<2/3<8/9<r. 
The above is symmetric in x and y, which gives another three terms. 

We are then left with four terms. For these we will use that, by the Mean Value 
Theorem, {xj{x) — Xi{y)Y ^ ll^XiliLk ~ Note that for the four remaining 
terms, 

X+Lt\+ , X-L[^\- , X+4"^X+ , X-i^"^X- , (5.4) 

we only need to take the supremum of |Vxj(C)l over the ^'s between x and y in the 
support of the relevant term. In this way we get: 

l(/,x±4"^x±/)l< // \f{x)\x±{^)\f{y)\x±{y)L'-r\^,y)d'xd^y 



< ^"jr ' / l/(^)lx±(x)((l/lx±)*G„)(x)d3x 

where Ga{x) = K2{a~^\x\) and cf(a) = sup|^|^2i I^Xj(2;)P- By first using the 
Cauchy- Schwartz inequality, then Young's inequality, we get 

\{f,X±Lfx±f)\ < ^^^^Il/X±l|2 II (I/I X±) * G„||2 < §^II/X±ll2 llG.lli. 
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Since ^ 

|jG„||i = J K2ia-^\x\)d^x ^ 4:7: ah^K2it)a dt ^ Gn^a^ 

(see l|A.6ll in Appendix ^ for t^K2{t)dt) we get the folfowing inequality: 

l(/,X±if X±/)l < ^Y^ll/X±ll^ (5.5) 

For two of the terms in (|5.4I) . x+l'^\+ and x+^a'^'^X+i this is sufhcient, since (see 
(Ol) : recall that I = a"^) 

c+ (a) = sup I Vxjf = c+a-2* , j = 2,3, 

\x\>2l 

and since t < 2/3 we get, using (|5.5|) . that 



5](m„x+4"^X+m.)<A^^c+«i-2* = o(a-4/3) , a ^ 0, 

as A'^ = \Z = XSa^^ (A and S fixed) and ||TOi||2 = 1- Similarly for x+L^'^Xa 
For the other two terms in l|5.4(l . note that 



ll/X-ll^= J \f{xr\x-i^)\'d^x^ J \ f{x)\'x-{^)d'x^{f,x-f) 

= {f,X-{xl+xl).f) = {xif,X-Xi.f) + {X2f,X-X2f) , 

since X- — X- ^^'^ Xi +X2 ~ 1 suppx-- Using this and H5.5|l . we obtain (since 
X-Xi = Xi): 

N 

53(m,,x-(Ll")+4"^)X-mO (5.6) 



Af A 



<Ca^ ^''(^^{Xlm^,Xlm^) + ^{x2m,,X-X2mi)^ 



where 

3 

C = :^(ci +C2) , cj-a-^- ^ sup |Vxj(x)|2 , j = 1,2. 

^ \x\<2l 

The two terms in (jS.fill will be estimated in the following two sections, the first 
one along with the energy at the nucleus, the second one with the energy in the 
intermediary zone. 

6. The energy near the nucleus 
In this section we estimate the energy at the nucleus, that is (see (I4.3f) '). the term 

N 

^(to„Xi{T(p) - V{x)+ap*(f,a * \x\'^}xim,). (6.1) 

i=l 

Also, half of the remaining term H5.6|l of the localisation error, treated in the previ- 
ous section, will be estimated here. We start by noting that p*4'a * is positive, 
so that we get a lower bound to 1)6. l|l by dropping this term. The remaining ex- 
pression will be treated by using the following result by Lieb and Yau 16, Theorem 
11] on the hyper-relativistic operator |p|: 
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Theorem 6.1. Let Co > and R > and let 

Hc,R = \p\~-\xr' -Co/R 

be defined on L'^(M.^) as a quadratic form. Let < j < q be a density matrix 
(that is, any bounded operator on L'^{M!^) which satisfies the operator inequality 
< 7 < q and for which Tr(7) < oo) and let x be any function with support in 
Br = {x \ \x\ < R}. Then 

TiixjxHcoR) > -4.4827 C4i?-ig{(3/4^i?3) J |^(^)|2 ^3^|^ 3) 

Note, that when x = 1 m Bji, then the factor in braces {} in (|6.2I) is L 

Here, TT{jh) is shorthand for '^)^{fk,hfk)^k, where {fk,lk) are the eigenfunc- 
tions and eigenvalues of 7. For more details, see Lieb ^^l- In our situation, q = 2. 
For our purpose, let 11 be the projection on spanjmi | i = 1, . . . N}, then 11 is a 
density matrix as above, and 



N 

' 1™« 

1=1 

Since suppxi C B{0, (1 + with 8/9 < r < 1, set = (1 + I3)a' and Co 

2(1 + /3)a'-i. Then 

5 



Tr(xinxi-ffco-R.) = '^{'>nt,xiHcQR.Xi'> 



T{p) - V{x) = + a-2 _ Q,-i _ 



2 

> IpI - a^^ |a;|"^ = HcoR + a~ 



since ^Jp^ + a ^ — a ^ > |p| — a ^ and 5 <2/tt. Including the first term in (|5.fcif) 
we now have, applying Ht).2|l . 



N N 



^(m,,xi{r(p) - V{x)]ximi) - C a" ^(wi, XiXi™z) 

i=l i=l 

N 

> ^{mi,Xi{HcoR + oT^ - C a^~^''}xiWi) 

i=l 
N 

> ^("^«, XiHcoRXim^) = Tr(xinxii?Coi?) > -Ca^--^. (6.3) 



The second inequality is valid for a small enough, since r < 1, so that a^*^^ ''^ ^ 
for a ^ 0. Since 3r - 4 > -4/3 (as 8/9 < r), the RHS of (gS is o{a-^/^), a ^ 0, 
which is the desired order. Note that the above procedure is what decides the scale 
a'",8/9<r<l, on which one can localise near the nucleus. 

7. The intermediary zone 
The energy in this area is given by the term (see (|4.3() ^ 

N 

^{'mi,X2{T{p) - V{x) +ap*(j)a* \xi\~^]x2mi). (7.1) 

i=l 
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The zone defined by the X2 was introduced to separate tlie outer zone defined by 
X3 and the support of the coherent states to be introduced later. As in the previous 
section we note that by dropping the term involving p * 0q * |a;|~^, we get a lower 
bound of the energy in (|7.1() . The remaining expression will be estimated by a 
generalisation of the Lieb-Thirring inequality (see Lieb and Thirring |14|'). proved 
by Daubechies in 01 page 518]. See also page 516 loc. cit. for the conditions on the 
function T{p). 



Proposition 7.1. LetF{s) ^ dt[T-^{t)f , whereT{p) = T(|p|) = ./\p\^ + cr 
a^^ as a function. Then 



N 

(V,^{T(p,)-F(a;0}^) > -qCj F{\V{x)\)d'x 



where C < 0.163. 



Note that this in particular means that the negative part of the spectrum of the 
operator T{p) — V{x) is discrete and that the sum of the negative eigenvalues of this 
operator is bounded from below by the quantity —qC J F{\V{x)\) d^x. To see this, 
let {ej}°^Q be these negative eigenvalues, eo < ei < . . . , and {gjjj^Q corresponding 
orthonormal eigenfunctions, and let tp be the Slater-determinant of the first N of 
the 5j 's. Then, by the above proposition, 

N 

~qC Fi\Vix)\)d'x < (V',^{r(p.)-^(2^.)}^> 

i=l 

N N 

= J2{g„ {Tip) - V{x)}g,) = ^e,. (7.2) 

Since the left-hand-side is independent of N, we get the statement by taking the 
limit TV — > oo. This will, as mentioned above, be used on the energy related to the 
cut-off X2, but also on the remaining half of the term x_(L^"'' + L^^^)x^ discussed 
in Sectional see (|5.6|l . First, let us calculate F: 



T{p) = T{\p\) - VbP + ^ T-\t) = ^¥+2(FH. 

Then 

a^N 3/2 



F{s)^ f {e + 2a-hf'^dt^ f {^f'\l^^) 
Jo Jo ^ 



dt. 



Now, by a Taylor expansion of the second term in the integral, we get 

(l + (a<)/2)^/^<l + ^i+^t2. (7.3) 

That is, for s > : 

ns)<{^''\\s^'^^'^/'^^p'^}. (7.4) 
The two terms we wish to estimate in this section are, as mentioned above, 

AT N 

'Y^{rni,X2{T{p) -V{x)]x2-mi) and Ca^"^'' ^(x2TOj, X-X2n^i)- 

i=l i=l 
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In order to do so, note that on suppx-X2 we have (x- being the characteristic 
function of B(0,2a'')) 

^ ' \x\- 20" - 

for a small enough, since r < 1, so that a^^'' —> as a ^ 0. Therefore, by the 
estimate (I7.4|l on F(s), and still for a small enough, we have 

N N 



^{mi,X2{T{p) - V{x)}x2mi) - C ^(m,, X2X-X2rn») 

1=1 1=1 

N 

>Y.^rn,,X2{T{p)~2V{x)]x2m{) (7.5) 



with V{x) = X2(a;)V^(a:^). Letting {ej,gj) be the negative eigenvalues and corre- 
sponding orthonormal eigenvectors for the operator T{p) — 2V{x) as before, we 
then have 

N N 

^{m^,X2{Tip) - 2V{x)}x2m,) >^{X2m,, ^^ej{g.j, •)5i}x2TO,) (7.6) 

i=l i=l j 

N N 

= Xm^jl(^2"ii,5j)p = Xl$I^jl("^*'^2gj)P >Y^ej\\x2gjf >J2^^- 

3 *=1 J i=l 3 3 

Here we used Bessel's inequality (remember that the m^'s are orthonormal), that 
e-j < and that < X2 < 1- Using on T{p) - 2V{x), in the limit N — > oo, we 
now reach (using H7.4II . (|7.5() . and H7.6|l ) 

iV TV 

X]('^»'X2{r(p) - V{x)}x2m) -Ca^~^''^{mi,X2X-^X2m) 

i=l i=l 

> -2C J F{2\Vix)\)d\ 

supp X2 

■2\3/2r2/^„^, ^,^5/2, 3a. ,,.7/2, 



> -c 

supp X2 

= -CAtt 
= ~C 



(^)"{|(2in^)i)^/^+^(2in.)i)^/^+^(2in.)i)^/^}^^. 



/2\3/2r2 2(5,5/2 , 3a. 2^,7/2 , a2 2,5,9/2-1, ,2 , 

4, i:z3 6(5, -('■+1) "(t + l) , 4(5^ l-3r l-3t . 

-(a2— Q;2j-|-_(Q- 2 —Q- 2 ) -|__^(q, 2 _ Q, 2 ) 



Since 8/9 < r < 1 and 1/3 < t < 2/3 , all of these terms are o(a~'*/^), which is the 
desired order. We note that it is this analysis that decides the scale a* of the outer 
cut-off X3- 
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8. The outer zone and Thomas-Fermi teory 
Up to order o{ar^l'^') we are now left with 

N 
i=l 

a rrp(^^3,,3^_^ 

This expression wih now be related to the semi-classical Thomas-Fermi energy. 
This is done by introducing coherent states, following Lieb and Yau in |15l proof 
of Lemma B.3]. Let g be the function chosen in Sectional that is, g S C(f'(M''), 
spherically symmetric, non-negative, supported in the unit ball -8(0, 1) of R'^ and 
such that / gix)"^ (Px = 1. Let ga{x) = a~^''^'^g{x/a''), 1/3 < t < s < 2/3, that is, 
4>a{x) = ga{xf' with a — a*. In this way, since N = XZ = X5a~^: 

which is also o(a~^/'^), a ^ 0. Define now the coherent states g^'*, p,q ^ by 

With T(p) the function ^Jp^ + — a~^, we then have the formulae 

a/) = (^// rf'pd'g(/,5r)(5r,/), 

{f:{V*gl)f) = jj d^d\V{q){f,gl:''m\f), (8.1) 

{f,T{p)f) > J I d^d\fip)if,gP'^){gP''iJ) - o(a-i/3). 

The proof of these formulae is carried out in Appendix ^ In this way, letting 
V{x) = S/\x\ — a p* \x\~^ (remember that (fias — g^): 

N 

^{■mt,X3{T{p) - V{x) +ap*(j)a* \x\^^}x3m^) 

i=l 

N ~ S 5 

= ^imz,X3{T{p) - V{x) * + |-| * - -TT Ixs^i) 

N 

= ^{mi,X3{T{p) - V{x) * tpa' ]x3m^) 

N 

^Wfll ^'^"^'9(^(?')-^(9))(El("^^X3,5r)l')-^o(a-V3). 

The second equality follows from Newton's theorem (since is spherically sym- 
metric): — *0qs = outside supp 0qs, and since supp X3risupp — 
for a sufficiently small (as s > t), 

IlK'X3{^*0a= - ^}X3"l») =0. 
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This is one of the reasons for introducing the intermediary zone by the function X2 ■ 
Note also that N o{a~^^^) — o(a~^/'^). Now, for a small enough, a**"* < 1/4, since 
s > i, so that if \q\ < ja* , then 



\x - q\ < \x\ < -a 



and so (wiXsj^Q^) since supp ga C 5(0, a^) and supp X3 C \ B{0, ^a*) 
That is, for a small enough 

suppJ(m,X3,5r)|' C M3 \ 5(0, 

so that for any ^ > we have, with M{p,q) — J2iLi l("^iX3i 5S''')P ^^'^ = 
max{±/, 0}: 

N 

^ jj d?pd\{f{j>)-V{q)){Y.\{m,x,,9l:')\^) 

1 ^ 



(27r)3 



> 



1 



(27r)3 



Tip) - (V(<Z) - 



\q\>iat 

since < M{p,q) < 1 and (xsm^jXa^i) < 11^^112 = 1- The first is seen by Bessel's 
inequality, since the m^'s are orthonormal and HxsSS'^lb < llsS'^IU = 1- In this 
way we have shown that for 11 > 0, p : R-^ ^ R and ^ e Hp = L'^{R^,C'^) : 

^ d'pd\[f[p)-{V(q)-ay) _ 



(27r)3 



(8.2) 

\x-v\ 

Choose now p to be the Thomas-Fermi density pj,p , that is, the function that 
minimises the Thomas-Fermi functional (here, 7 — (37r^)^/3); 



£tf{p) = \lj p{xf" d\ ~ J p(x)^^ + ^ // 



P{x)p{y) ,3 ,3 /o 

-; T-dxdy (8.3) 



{peL^/3(M3)nLi(R3)|p>0, / p{x)d\<N^ 



over the set 



(For the Thomas-Fermi theory, see Lieb and Simon [13] and Lieb Then p^/ 

satisfies the Thomas-Fermi equation; 

Z . 



2/3 



for some unique p = p{N). Furthermore, 



(8.4) 



for < Z: 
for N > Z: 



p^p{x) d^x = and p{N) > 0, 



p^/(a;) d'^x = Z and A*(^) = 
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(see Lieb and Simon Theorems 11.17, 18 and 20]). In this way, / pj.p {x) (Px < 
N imphes N > Z, and therefore /i(A^) = 0, so that we always have 



fi{N) I p^'p{x)Sx^^i{N)N. 



(8.5) 



Let £tf{N, Z) = £tf{Ptf) ^^^^ define the Thomas-Fermi potential by 



V^/ix) EE Z/\x\ - p^>" * \x\-' - piN), 

then we have the following scaling f [Tgl (2.24) p. 608]) (remember, that A = N/Z is 
fixed) : 



£tf{N,Z) = Z'^^£tf{X, 1) = -CTF(A)Z^/^ 
V^/ix) - Z^/'V^p\z'/'x) = Z^/'Vtf{Z'/'x). 



(8.6) 
(8.7) 



The idea is now to estimate the difference between the integral in (|8.2() (with 



P = Ptf ^^'^ P ~ ^'^d 
a 



(2n) 



d'pd\[^-{^^-p^F'*\<ir-piN)) 



\i\>ia* 



This is done in two steps: first, we change the domain of the integration, then we 
change the integrand, each time estimating the error. 
First, 



1 



(27r)3 



|g|>ia' 



d'pd\{f[p)-aV^f{q)) + 



1 



(27r)3 



d^pd\{T{p)~aV^Y{q)) 



(27r)3 

Since t(p) > 0, we get 



1 



(27r)3 

|«|>ia'; fipXaV^ji.'' {q)<a 



d'pd\{T{p)-aV^Y{q)). 



d'pd\{aV^/{q)~Tip)) <a 



d^pd\V^f{q). 



|g|>iQ*; f(p)<aV^f(q)<aS± 



Using the scaling (|8.7(l and the change of variables to — S'^/^a ^^^g, the above is 
equal to 



51/3^,2/3 



dyd'^uoVrFi'^)- 



|^|>l5l/3„t-l/3 

f{p)<S^''^a-^''^VTF{i^)<a4- 

The limits in the integral means that 

2^4/3^-4/3^^^ < p2 < 25^/^a-^"'VTF[uj){l + ^S^'^a^'^'VTFi^)) 
so that with 

X = 25^'^a~^/^VTF[u;) , Y - ^S^/'a^/^VTF{uj) , Z ^ \p\^ , W = ^S^/^a'-^/^ 



16 



THOMAS 0STERGAARD S0RENSEN 



we have 



^ ^ (4^)2ji/3a2/3 / d\u:\ \oj\''Vtf{uj){ / ^ dZ 



w 

oo 



^3/2 



= (4^)2<5V3„2/3 / d\u;\\i0\^VTF{^)^{{l + Yf/^-l) 



w 



By the Taylor-expansion lf73|) . we have (1 + F)^/^ < 1 + |y + fF^, and so 
Jw 



Usmg that [x) < Z/\x\, smce /i(A^) > and pi^p > (remember that Vtf 



N,Z 



V^p)^ we arrive at 



w 



-5/2 



W 



a"2/3^-l/2 _^ ^-3/2 ^ ^-2/3^1/6-t/2 ^ ^(l-3t)/2 

o(a-5/6) + o(a-i/2) 



since i < 2/3. This means, that 
1 



(27r) 



> 



(27r)3 



d^pd^g (f (p) - aV^f{q)) - o{a-^/') 



\<i\>i 



'-<aVj,ji (?) 



Next note that since \q\ > ja* and aV^j^ (q) < S/\q\ in the area of integration, 
we here have that 



? 4 
P „3P 



z o 



In this way, we get 
1 



(27r)3 



1 



> 



(27r 



l3 



d^pd\{a^-a'^-aV^'{q)) 



1 



(27r)3 



^d^pd\. (8.9) 
8 



|g|>ia*; aii-<ay~/(9) 



|g|>ia*; aii-<aV~^'^(g) 
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Note that 



II d^d\{at^-aV^f{<l)) 



(27r) = 



> 



(27r)= 



Let us now look at the last term in 1)8.9(1 . Again using that {x) < Z/\x\, we 
have that 

^d'pd\ < J J ^Spd\ 

^ypIzM |^|4 



(47r)2 / d\q\ \q\' 



27r2(2Z)7/2 



7 Jia' ' 



-t/2 



Using this, we then get the following 
1 



(27r)3 



d^vd\{f{v)-oiV^f{^)) 



l9l> 



%-<QyyJ, (g) 



> -■ 



(27r)3 



(2Z)^ 



Hence, since (5 = is fixed and i < 2/3, we have 



a ^ 0. 



Summing up, we have now proved that for ij; ^Tip = /\^ L^(R'^, C^] 



(27r)3 

a 
2 



d^p d\ 



a ff Ptf {x)Ptf (y) ^3„j3, 



\x - v\ 



d^xd^y-afi{N)N -o(a-^/^). (8.10) 
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Integrating firstly in p in the first integral in (|8.1()|l . we get, for each q fixed: 

"2 Z 



2 (j^\-Ptf'\<i\-'-KN)) 



16%/27r 



15 



(8.11) 

5/2 



O ^ Vrj- 



The [•••] + , since, if the term in brackets is negative, the integrand in (|8.11|l will 
be zero. 

Now, because p^p satifies the equation 18.4|l . we get, that 



[ Z 

7\ 



N,Z I 1-1 

Ptf * \l\ 



5/2 



3/2 N.Z/ ^ 

l' Ptf W 



191 



In the last equation, no [■•■]+ is needed, since, if the last term is negative, p^p is 
zero, because of (|8.4(l . In this way, by the above and by H8.5|l : 

.2 7 



(27r)3 



d^p d\ 



2 



(t 



a Ptf {^)Ptf (V) .3^-3, 



N,Z I 1-1 

Ptf * l^l 

N.Z, 



\x - y\ 



d^xd^y ~ a^i{N)N 



. N.Z, 

"^7 / Ptf 



+ « / p'Mq) Ptf * W d'q + apiN) I p'^p^{q) d'q 



N.Z / 



N.Z / \ i\ ,Zj ( \ 

Ptf \^)Ptf KV) ,3^ ,3, 



N.Z I 



V - y\ 



d^'xd^'y - ap{N)N 



1 P^F'{^f''d\- p^p^[x)^d\ 



N.Z, \ N.Z, \ 

Ptf \^)Ptf \y) 
\x - y\ 



d^x d^y 



\x\ 2 
= a£TF{N, Z). 

Since Hrei = a^^H, and Z = 6a^^, with S fixed, < S < 2/7r, this shows, that for 
aU^e = A^i'(K^C2): 

{iP, Hrei^) > - CtfZ'"!^ - o{Z''l^) , Z ^ 00, 

because of the scaling (|8.t)|) . This ends the proof of Theorem 11.11 □ 



Appendix A. A formula for the kinetic energy 

In this appendix we shall prove the localisation-formula for the operator 
iyp^"+cr^ (which is the equivalent of the IMS Localisation Formula for the Laplace 
operator —A, see Cycon, Froese, Kirsch and Simon J?, Theorem 3.2]). Let firstly 
Ki be a modified Bessel-function of second order, defined on (0, 00) by 
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It is easily seen that K2 is well-defined, decreasing and differentiable. Other prop- 
erties of K2 will be derived later. Let then Xj, j = I, . . . , fc be smooth positive 
functions on M'^, such that Xjix) = 1 for all x in R-^ and define on L^(M'^) the 
bounded operator L*^"^ by the kernel 

Then for / G one has the formula: 

k 



if, Vp' + a-\f) = ^(/,X,v^^+^X,/) - (A.l) 

The proof of the localisation formula (jA.ip will be a consequence of the following 
formula: 

Lemma A.l. For f e S{R^) , 

(A.2) 

Proof. Let / be the Fourier transform of /. Note that by dominated convergence 
in momentum space, we have 

(/, ^p^ + a-^f) = liin ^{(/, /) - (/, e-Vp^+"-V)}. 



To calculate the integral kernel exp[— t-^/p^ + a ^](a;, y), expand the Fourier trans- 
forms: 



This is justified by the fact that / £ 5(]R3). Now, for x, y fixed, choose polar 
coordinates 6, 0), for p such that {x — y) ■ p — — \p\\x — y\ cos9. Then 

00 p27T PTT 



"'0 "'0 



g-Vp'+a ^^-t\p\ \x-y\ cos S q j^p ^j^j 

27r^ \p\^e'*^P'+"''(^J e'IPll^-2/l«dw^ , u = - C( 
/ |p|e-Vp^+-%in(|p| |x - y|) d\p\ 

F ~ yi Jo 

-^ta-^\x - y\{\x - + t^)" - yp + i^)!/^] 
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where the last equahty is given in Erdelyi, Magnus, Oberhettinger and Tricomi 
p. 75, 2.4 (35)]. In this way, 



Now, letting Ft{p) = e *Vp^+" the above shows that 

1 /■ . , /2" , X2[a-i(|a;P+t2)i/21 

Ft{x)^-—rj. Ft{p)e''--PSp=J-ta-^ n VM ; J 



(27r)3/2 7 '''^^ V TT |a;|2 +t2 

and therefore, for all y e R'^: 

Hence we get, using (jA.3|) and (|A.4|I . which are both symmetric in x and y, that 

^{(/,/)-(/,e-Vp^+""V)} 

= ^{(/,/) - + |{(/,e-*""V) - (/,e-Vp^+"-V)} 

(/,/) 



i- 



/ ^ ((1/(^)1' + - - 



to-2 j^2[«-^(|x-;/p+t2)i/2] 

Cancelling t and noting that 

hm = ^(e-*" ) 



= —a ^, 



we get that 



limi{(/,/)-(/,e-'Vp"+<.-V)} 



This proves the lemma. □ 
Now, to prove the formula l|A.l|l . we simply use the fact that x^ji^) = 1 for 



^|x,(a;)/(x)-x,(y)/(y)l^ 



|/(x) - /(y) |2 + ^ (x) (7M/(^) + W)f{y)) (x, (^) - X, (y)) . 
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Note that Xj/ G '5(K'^), since Xj is smooth and bounded, so that using the formula 

k k 



, -2 



(27r)2 



J \x y\ 



Using now that 



K2{a ^\x-y\ 



X, {x)f{y)f{x) {x, {x) - X, {y)) 'V _' I, d\ d'y 

p y\ 

K2{a-'\x-y\) 



\x y\ 

simply by interchanging x and y, we finally get from l|A.5|) that 



{2nrJJ "-'^^-^^ ' -^-'^ \x-y\^^ 

which, using l|A.2p . proves the formula (|A.ip . □ 
We now derive two facts about the function K2: 

1"°° "^TT 

t^K2{t)dt = ^, (A.6) 
i^2(i) <4yje-*(l + l + ^) foraU<eM+ (A.7) 
The proof of (|A.6|1 is straightforward by using the definition of K2'. 

POO -I />OC 

t^K2it)dt^ ei- xe-i'^''+''''Ux) dt 

= LJ ^(^J ee-'-^^+^'Ut^dx 

where the interchanging of the order of integration is allowed by Tonelli's theorem. 
By applying partial integration three times, 

[x + x ^) 
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and so 



1 f"" 16x f°° 
t^K2it)dt^-l :r-dx^4: -5- da; 



4 ^ L _ ^ 

(z2 + i)3'7 """32^ ~ 2 ■ 



4 ■ 27rz Resf , „ — nr-M = Stti 



For the estimate (|A.7|I . we need to rewrite if 2. This is done foUowing Gray and 
Mathews |7| pp. 50]. 



Observation A. 2. 

m) = xliT^'!-' I e-ie'Hl + hr^K. (A, 



To prove the observation, we start on the right-hand-side of IjA.Sp . Setting 



t-\- ^ = ^t"^ + T], one gets, since then 77 = _|_ 2t^^ that 



RHS 



1 /"°°g-^/^^/'^^3/2_d?7_ 



Using the formula 



^ 2a 



(which holds since both sides satisfy the differential equation df /db = —2a/, f{b 
0) = y/Tr/2a) with a = ^/t^ + r], b = 1/2, we arrive at 

g-(t^€^ + l/(24) = )^-5^^ 



(2t)2 

since one has the formula 



poo 

/ e-^'^\'/' dr, ^ r'ni) 

Jo 



Making the change of variables x = , we finally get 



RHS 



1 f 

2 Jo 
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Now, to prove the estimate (|A.7|I , use the Tayloer expansion ^ on the integrand 
in (|A.8|I . to get 



2 

2t r(|) ^ ' ^ ^ 



3 

it 



POO n poc 

(r(l) + |r,i, + 5|,r(|)) 



2i rxl) 



TT 15 105 , 1 1 , 

2^ ^ (1 + ¥ + 128^) ^ ' + 2^ + 

Appendix B. Introducing coherent states 

In this section we will introduce coherent states and prove the formulae in section 
O The error introduced by using coherent states will also be estimated here. 

Lemma B.l. Let g G C5^(M^) he spherically symmetric, non-negative, supported 
in the unit ball and such that \\g\\2 ~ 1, and let = g{x — q)e^P^ . Then 

(/,/) = (^// d'pd\{f,g'^-'^){gP^''J) 

{.f,{V*\g\')f)^j^ jj d'pd\V{q){f,gP^''W-\f) 
1 



(/, Vp^T^/) > ^ j J d^d\^p^ + a-^f,gP'^)igP''>J) 

-3a||Vg||LVol(supp g)\\f\\l (B.l) 

Proof. The idea of the above formulae is to write the identity and other operators 
on i^(]R^) as superpositions of the one-rank operators Wpg = ( , g^''^)g^'^. To prove 
the above formulae, start with the right-hand-side of the second formula (the proof 
of the first formula is similar, just more simple): 

^ jj d^pd\V{q){f,g^''^W^M) (B.2) 

1 



d^pd\V{q) / f{y)g{v ~ q)e-^Py dhj / f{x)g{x - q)e-'P- d\ 



(27r)3 

Notice, that the function in the last brackets is (27r)"^/^ times the Fourier-transform 
of the function Fq{x) = f{x)g{x — q). In this way we get, by Parseval's formula: 

^=11 d^pd\V{q)\F,{p)\' = J d\V{q)\\F,\\l^ J d\V{q)\\F,\\l 
d\V{q)[ I \f{x)\'\g{x-q)\'d'x) 

'd^x\f{x)\'[ I V{q) \g{x-q)\'d\)^{f,{V*\g\')f). 
This proves the second (and the first) formula. 
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To prove the formula for the operator ^Jp^ + a ^ , note that 
j g{x - qf d\ = 1 for all x in M^, 



so that, by the symmetry of the operator yp^ 



if, y/p^ + a-^f) = - 1 1 fix)g{x - qfiVp^+^f) (x) d\ d\ 



{^p^ + a-^f){x) g{x - qff(x) d\d^x 
\jj W)9,{xf{^p^ + a-y) (x) d\ dh 
+ J{x)(^y/p^ + a-^{g,'f)){x)d\d'x. (B.3) 



Here, gq{x) = g{x — q). Remembering that gq{xY is reel and letting gq^ denote the 
multiplication operator defined by this function, we have 



131= {gqf) {x) [ vV+^(5g/)] {x) d\ d\ (B.4) 



+ 2 / / /(■^) (59 Vp' + + ^P^ + «-'.9,' - + a-^gq)f {x)d\d\ 

__ I 
" 2 



f{x){Lqf){x) d\dh 

+ JJ {J V^+^{ J e-^Pygq{y)f{y)d\)e^P^d^)gq{x)J{x), d\d\ 
where 

{Lqf){x)= (B.5) 
/ { / [9,{yf + 9,{xf - 2gq{x)gq(y)\ v^^T^e^^^^"^) dPp]f{y) dV 
The second term in (|B.4|I is equal to 

J J d^d\^Jp^T^(^jl(x)9,{x)e^''^d\^ (^j f{v)gqiy)e-'Py d^y^ 



d^d\V^T^{f,g''''>W'^f). 



The first term in (jB.4(l is the error, which will now be estimated. Keeping x and y 
fixed, we have, as showed in the proof of (|A.2|I : 



Lqix,y)^ j [gqiyr+gqixr-2gqix)gqiy)]Vp^T^e^P^^-yU'p 
= [9q{x) ~ gq{y)] 

In this way, using the same ideas as in Sectional we reach the estimate 

Lq{x,y) < llVg^ll^ ^^2(0 ^\x - y\) (xsuppg, (a;) + Xsuppg, (y)) , 
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where Xsuppg, is the characteristic function of suppf;^. This gives us that 
jLq{x,y) d\< j llVg^ll^ ^ii:2(a"V- y|) (Xsupp5,(a;) +xsupp.gg(y)) rf^g 

= 2 ||V.9||2^ ^ K2{a-^\x - y\) Vol(suppg). 

By this we finally get, by using first Cauchy- Schwartz's, then Young's inequality, 
that 



< 



f{x) j Lg{x,y)d\f{y)d^xd^y 

_2 

\fix)\ (2 \\Vg\\l ^K2{a-'\x- Vol(supp <?)) \f{y)\d\d\ 



< 2 IIV5IIL 4^^11/112 II I/I *G„||2Vol(supp5) , G^{x) ^ K2{a~'\x\ 



^ l|V.g||^^ 11/11^ ||G„||iVol(supp5) 

— 2 

= l|Vg||L^67rV||/||2Vol(suppg) (seeElfor ||G„||i) 
= 3a||Vg||^Vol(suppg)||/||i 



□ 



For the case (|8.1ll in Section |S| let the coherent state 5^'* be defined from the 
scaled version of the function g chosen there — that is, g € G^(R'^), spherically 
symmetric, non-negative and with support in the unit ball -6(0, 1) of M'^. Then the 
coherent states are 



ipx 



In this way, ||V.ga||^ a ^^IIV^H^ and Vo1(supp5q) = ^a^^, and therefore 

since, as s < 2/3, 

3aa-5^ llVgllLf 11/11^ = C a'-'^ = o{a~'/') , a - 0. 
This proves the formula H8.1|) . since 

(/,/) = jj d^d\{f,gP^'^W^\f) 

and T(p) = ^Jp^ + a^'^ - oT^ . 
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